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ABSTRACT 

The nucleosynthesis and ejection of radioactive 26 Al (ti/2 ~ 0.72 Myr) and 60 Fe, (tj/2 ~ 2.5 Myr) 
into the interstellar medium is dominated by the stellar winds of massive stars and supernova type II 
explosions. Studies of meteorites and their components indicate that the initial abundances of these 
short-lived radionuclides in the solar protoplanetary disk were higher than the background levels 
of the galaxy inferred from 7-ray astronomy and models of the galactic chemical evolution. This 
observation has been used to argue for a late-stage addition of stellar debris to the Solar System's 
parental molecular cloud or, alternatively, the solar protoplanetary disk, thereby requiring a special 
scenario for the formation of our Solar System. Here, we use supercomputers to model — from first 
principles — the production, transport and admixing of freshly synthesized 26 Al and 60 Fe in star- 
forming regions within giant molecular clouds. Under typical star-formation conditions, the levels of 
26 Al in most star-forming regions are comparable to that deduced from meteorites, suggesting that 
the presence of short-lived radionuclides in the early Solar System is a generic feature of the chemical 
evolution of giant molecular clouds. The 60 Fe/ 26 Al yield ratio of k, 0.2 calculated from our simulations 
is consistent with the galactic value of 0.15 ± 0.06 inferred from 7-ray astronomy but is significantly 
higher than most current Solar System measurements indicate. We suggest that estimates based on 
differentiated meteorites and some chondritic components may not be representative of the initial 60 Fe 
abundance of the bulk Solar System. 

Subject headings: ISM: kinematics and dynamics — stars: formation — nuclear reactions, nucleosyn- 
thesis, abundances — protoplanetary disks — meteorites, meteors, meteoroids 



1. INTRODUCTION 

Due to the low opacity of the Milky Way to 7-ray 
emissions from the radioactive 26 Al and 60 Fe nuclei, 
7-ray observations provide estimates of the mean 60 Fe 
and 26 Al abundances in the Galaxy. Combined with 
galactic chemical evolution models, these observations 
can be used to determine the steady-state abundance 
of 26 Al and 60 Fe in the interstellar medium and, im- 
portantly, a robust estimate of the 60 Fe/ 26 Al ratio of 
the Galaxy (jDiehl et al.ll2006t I Wang "et^ [2001 . Mete- 
orites and their components, including the Solar System's 
first solids — calcium-aluminum-rich inclusions (CAIs) — 
contain evidence for an early presence of 26 Al and 60 Fe 
at levels that were apparently higher than expected from 
the galactic back ground abu ndance, especially for 26 Al 
(jKrot et al.H2009t IHuss et al.ll2009f ). Moreover, the vari- 
able initial abundances of 10 Be, a short-lived radionuclide 
uniquely synthesized by spallation reactions, in CAIs 
recording uniform initial 26 A1/ 27 A1 values indicate that 
no more than ~2.5% of the solar system's 26 Al inven- 
tory m ay have been produced by local irradiation pro- 
cesses (|Kennicutt fc Ev ans 2012). These observations 
are interpreted as reflecting a late-stage contamination 
of the nascent S olar System from a nearby supernova 
(jGoswamil 12004( 1. However, the majority of Sun- like 
stars form in cold and dense molecular cloud regions 
encap sulated in Giant Molecul ar Cloud (GMC) struc- 
tures (|Kennicutt fc Evans! 120121 ) . an astrophysical envi- 
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ronmcnt different from that of the average galactic in- 
terstellar medium (ISM). Therefore, an overall galactic 
abundance level of 26 Al and 60 Fe is not necessarily rep- 
resentative of levels encountered in star-forming regions, 
where the cold and dense gas out of which stars form may 
inherit material enriched in freshly synthesized matter by 
nearby, recent generations of stars. 

To better understand the abundance levels of 26 Al and 
60 Fe in star-forming regions, we used three-dimensional 
magnetohydrodynamic models of self-gravitating GMCs 
to simulate the time-integrated production and ejec- 
tion of these radioactive species via the supernova and 
stellar wind mechanisms of ensembles of stars, track- 
ing their incorporation into star-forming clumps. In 
these simulations, the ISM is assumed to be driven by 
the kinetic energy originating from supernova input, as 
well as from larger, galactic-scale, density perturbations, 
producing giant molecular clouds with morphology and 
kinematics controlled by gravo-turbulent fragmentation. 
A hierarchical structure results, with GMCs contain- 
ing smaller scale cold Molecular Cloud (MC) fragments, 
which in turn contain filamentary structures where suf- 
ficiently dense and massive pre-stellar cores may be- 
come gravitationally unstable and collapse into new-born 
stars, with initial masses di stributed according to an 
Initial Mass F unction (IMF. iPadoan fe Nordlundl 120021: 
Chabrier 2005). Sufficiently massive stars explode as su- 
pernovae after a few million years, enriching their sur- 
roundings with freshly synthesized metals as well as con- 
tributing kinetic energy back into the ISM. 

2. METHODS 

We use periodic simulation boxes with sizes rang- 
ing from 40-80 pc, filled with self-gravitating, magne- 
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tized gas of average mas s density con sistent with Lar- 
son's mass-size relation (|Larsonl 119811) : this results in 
H2 number densities between 30 and 15 cm -3 , cor- 
responding to total masses ~ 10 5 and 4 ■ 10 5 M Q 
respectively. Experiments are initialized with the 
unigrid magneto-h ydrodynamics STAGGER code (see 
IKritsuk et al.ll20lTI ) that includes statistical star forma- 
tion and individual supernova ex plosions, with tracking 
of yields (jLimongi fe Ch icffi 2006) using a passive scalar 
transport equation. The passive scalar variable is as- 
signed the decay time of its representing isotope ( 26 A1 de- 
cays to 26 Mg with a half-life of 0.72 M yr and 60 Fe decays 
to 60 Ni with a ha lf-life of 2.62 Myr; iNorris et al.lU983l; 
iRugel et al.l I2009D . Solar abundances ( Asplun d et al.l 
l2009h are assumed for the stable reference isotopes 56 Fe 
and 27 A1. Bit-wise identical simulation runs are executed 
in pairs, differing only in half-life and yields. Because the 
passive scalar does not affect the dynamic evolution of 
the simulation short-lived radioisotope (SLR) data from 
different runs can be superimposed. 

The periodic boundary conditions represent a GMC 
fragment surrounded by similar or larger amounts of 
nearby gas, with the surrounding gas pressure prevent- 
ing break-up due to supernova explosions and other feed- 
back processes; GMC life times are consiste nt with the 
'star f ormat ion in a crossing time' paradigm (|Elmegreenl 
2000, l20ll . The view that feedback is very impor- 
tant for cloud breakup, while popular, is actually not 
suppor ted by convincing, l arge scale numerical experi- 
ments (IKritsuk fe Normanl 12011). Inertial forcing from 
larger scales is replaced our mode l by solenoidal forc- 
ing in a shell of low wavenumbers ()Padoan fe Nordlundl 
12002ft . here 1/L < k/2n < 2/L, producing velocity dis- 
persions of the cold and d ense medium consistent with 
Larson's velocity relation ([Larson! [T981). Heating and 
cooling is modeled as schematic heating by UV-photons 
(|Osterbrockl I1989T) quenched i n den se gas according to 
the recipe of iFranco fe Coxl ()1986l ). and an optically 
thin c ooling function consistent with iGnedin fe Hollonl 
(2012). Connected cold and dense regions of sufficient 
mass (> 10 3 Mq for occasional formation of 10 2 Mq 
stars) are converted to stars at a given star formation 
rate, g, with an IMF-dist r ibutio n of masses consistent 
with iPadoan fe Nordlundl ()2002h . The star formation 
rate per unit free fall time e was varied from 0.02 to 0.1, 
covering a range con sistent with observations, numerical 
models, and theory dMurravl 120111 : IPadoan fe Nordlundl 
l201ll: IkTWioIz fe McKedl2005H . 

Stars > 8 Mq are followed until they explode as 
Type II supernovae , after lifetimes interpolated from 
iMevnet et al.l (|1994ft . Supernovae are assumed to release 
10 51 erg of thermal energy, distributed in a mass corre- 
sponding to the entire mass of the progenitor star. We 
combine the small expected integrated wind yields into 
the generally larger supernova yields. Supernovae explo- 
sions result in initial temperatures of ^10 8 K and initial 
expansion velocities of ^1000-3000 km s , which for 
brief periods of time enforce time steps ~years. To allow 
the initial state of constant density and zero velocity to 
develop into realistic GMC initial conditions, we let the 
forcing generate a supersonic veloci ty field with av erage 
velocity dispersion consistent with iLarsonl ([1981), and 
wait until density fluctuations generate a first generation 



of stars, where the most massive stars already begin to 
explode as supernovae. 

When supernova driving of the medium becomes sig- 
nificant a snapshot is passed to a locally-modifie d version 
of th e adaptive mesh refine ment code RAMSES (jTevssierl 
I2002t iFromang et al.l 12005) . with stellar lifetimes, su- 
pernova yields, and ISM heating and cooling consis- 
tent with that used in the STAGGER parts of the ex- 
periments. In RAMSES experiments there is no external 
driving and star formation is handled individually, with 
new-born stars represented by accreting sink p articles 
(jPadoan fc Nordlundl l20lTI : IPadoan et~al1 127)121 ). The 
RAMSES code allows for tracking the transport of sev- 
eral nuclides in a single simulation. Massive stars are 
tracked until they explode as Type II supernovae. We 
use a 128 3 root grid and refine to a local resolution cor- 
responding to 4096 3 (cell size at the highest refinement 
level As « 0.01 pc w 2000 AU). Tests with 256 3 RAMSES 
root grids and corresponding STAGGER code resolutions 
did not give significantly different results. 

3. ENHANCED ABUNDANCES OF 26 AL AND 60 FE IN 
EVOLVING GMCS 

We find that ~10 Myr of evolution is needed for 
the model GMC to develop properties similar to young 
GMCs, where most of the massive stars are waiting to 
explode as supernovae. Over the next ~20 Myr, the 
results show a gradual increase of supernova feedback, 
with levels of both kinetic energy and average yields be- 
having in a similar way. This time interval represents 
the lifetime of the model GMC where we track the pro- 
duction, transport and admixing of stellar-derived 26 Al 
and 60 Fe into cold and dense star-forming gas. Over 
this ti me period, con sistent with estimates of GMC life- 
times (jMurravl 120111 ). the average concentration of 26 Al 
and 60 Fe in star-forming gas increases, with 26 Al/ 2 Al 
and 60 Fe/ 56 Fe values ranging from ^10~ 6 to ^10~ 4 and 
^10~ 7 to ^10~ 5 , respectively (Fig. [1]). There is signif- 
icant variability in the 26 A1/ 27 A1 and 60 Fe/ 56 Fe ratios 
amongst the individual star-forming clumps interspersed 
within the GMC throughout its evolution. A first im- 
portant quantitative result is that the concentration of 
26 Al and 60 Fe in dense star-forming gas is systematically 
higher than the inferred galactic background abundances 
of these radioactive species, estimated to 26 A1/ 27 A1 and 
60 Fe/ 5 6 Fe values of ~ 2. 3 • 10~ 6 and - 1.6 • 10~ 8 , respec- 
tively (jHuss et all 120091) . Transient bursts of 26 Al and 
60 Fe abundances observed in the first 10 Myr of the sim- 
ulation predominantly reflect single supernova events of 
the most massive stars (M Q >20). In later stages of the 
simulation (>10 Myr), most of the 26 Al and 60 Fe budget 
is contributed from the longer-lived lower mass massive 
stars (8<M Q <20), typified by more frequent lower yields 
compared to their rarer and more massive counterparts. 
This results in a more flat distribution of 26 Al and 60 Fe 
amongst star-forming regions. 

Figure [2] shows RAMSES results for the 5-10 Myr in- 
terval of the model GMC, computed with a finest grid 
resolution of 0.01 pc, which allows us to trace the con- 
centration of 26 Al in cold and dense gas accreting to in- 
dividual stars. The average 26 Al abundance during that 
time period is comparable to that obtained with the low- 
resolution STAGGER code although the variability of the 
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Fig. 1. — Supernova mass as a function of time (top), and values 
of the ratios 26 A1/ 27 A1 (mid) and B0 Fe/ 56 Fe (bottom) in dense, 
star-forming gas of an evolving 40 pc, 10 Mq GMC with e = 0.05. 
Grey dots show values sampled in dense, star-forming clumps, with 
10% and 90% percentiles marked in blue an d cyan. Dashed lin e 
is the modeled bulk solar value of 2.8 ■ 10 -5 (Larscn ct al][20lj). 
Time is indicated in Myr after formation of the GMC (10 Myr 
offset from simulation start). 

26 A1/ 27 A1 values is greater in the high-resolution RAM- 
SES case, which is able to maintain larger amplitude in- 
homogeneities inside the compact molecular cloud struc- 
tures. We note that gas accreting onto our 'sink particle' 
representation of stars do so at densities ~ 10 7 cm" 3 . 
These contributions shows similar levels of 26 Al — within 
the rather large statistical spread-as that of the average 
star forming gas, sampled at ~ 10 4 cm~ 3 . This illustrates 
that, although inefficient, star- formation in cold molec- 
ular clouds is a locally rapid p rocess, occurring on time 
scales of order crossing times (jElmegreenl 120001 ) . Thus, 
we infer that the enhanced abundance of 26 Al and 60 Fe 
in star-forming gas observed in our simulations — as com- 
pared to the average galactic abundances — is a generic 
feature of the chemical evolution of GMCs. 

Figure [3] shows the ratio of 60 Fe to 26 Al in star-forming 
gas and individual supernova ejecta, as well as the ratio 
of accumulated yields at any one time. The latter pro- 




time [Myr] 

Fig. 2. — Similar results as those shown in the middle and bottom 
panels of Fig. [T] but computed with the RAMSES adaptive mesh 
refinement code. The timescale corresponds to that of Fig. [T] The 
slanted gray lines in background indicate the decay rates. Yellow 
thin curves show instantaneous relative 26 A1/ 27 A1 abundance in 
gas accreting to sink particles representing individual stars. 

vides a direct estimate of the ratio of expected 7-ray 
line fluxes, since the large majority of ejected radioac- 
tive 26 Al and 60 Fe nuclei will have had time to decay. 
For both the STAGGER and RAMSES simulations, the pre- 
dicted 60 Fe/ 26 Al ratio of total yields is slightly less than 
0.2, consist ent with the obser ved 7-ray line flux ratio of 
0.15±0.06 (|Wang et alJl2007t ). The instantaneous ratio 
of average number densities of 60 Fe and 26 Al is larger 
by the ratio of their half-life, 3.6. Moreover, the aver- 
age 60 Fe/ 26 Al ratio in star-forming gas is expected to 
be larger than the global average, reflecting the time 
required for hot supernova ejecta to be converted into 
dense star-forming gas. Collectively, this implies that 
the expected average 60 Fe/ 26 Al ratio in star forming gas 
is typically considerably larger than 0.2 (Figure [3]). 

4. TIMESCALES AND MODE OF ADMIXING SUPERNOVA 
EJECTA INTO STAR-FORMING GAS 

The current perception of transport and mixing of 
chemical yields in supernova-driven ISM is that this pro- 
cess is rather inefficient, leading to mixing timescales 
that are longer than the t ypical lifetime of GMCs 
(|de Avillez fc Mac Low! [2002D . However, these simula- 
tions are based on the passive seeding of particle tracers 
to a turbulent medium and, therefore, the tracer particles 
are unrelated to the supernova events driving the ISM. 
In contrast, our code allows for linking passive scalars 
and particle tracers used as proxies for 26 Al and 60 Fe 
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Fig. 3. — 60 Fe/ 26 Al ratios in the star-forming gas of an evolving 
GMC for cold star-forming clumps (grey) and individual super- 
nova yields (orange) , from low-resolution STAGGER code simulations 
(top) and from high-resolution RAMSES simulations (bottom), with 
the ratio of accumulated yields (red) . 

with the kinetic energy of individual supernovae, thereby 
providing realistic timescales for the transport and ad- 
mixing of supernova ejecta into star-forming molecular 
gas. For example, transfer of kinetic energy to the gen- 
eral ISM from slowed-down SN ejecta drives additional 
compression and conversion to turbulent, cold MC gas, 
which imparts collapse of some gas fraction within its 
dynamical timescale. To illustrate the rapid incorpora- 
tion of supernova-derived 26 Al and 60 Fe into star-forming 
regions we show in Figure [4] accumulated fractions of ini- 
tially hot supernova ejecta that reached cold, dense, and 
potentially star-forming gas, after a given time. These 
results indicate that in some cases (of the order 1%; cf. 
Fig- HI), only a few 100,000 years are needed for part of the 
supernova yields to be incorporated into self-gravitating 
cores. A significant fraction of supernova ejecta (5- 
10%) has been transformed into star-forming gas over 
a timescale of only 1-2 Myr. 

To better understand how stellar yields are transported 
into star- forming regions, we show in Figure [5] examples 
of projected total and 26 Al number densities during the 
expansion of a supernova bubble over a time interval of 
1 Myr. Star-forming molecular cloud fragments are visi- 
ble as filamentary structures in the total number density 
projections, and these do not change dramatically dur- 
ing the 1 Myr period. In contrast, the supernova bubble 
expands to cover a significant fraction of the gas volume. 
As illustrated by the right hand side panels, where the 
upper two panels have the same absolute normalization, 
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Fig. 4. — Accumulated histogram of the fraction of supernova 
ejecta trace particles (N) that has reached a cold and dense (N co id), 
potentially star-forming state after a given time. The solid curve 
is an average over the six individual cases shown by dotted lines. 
These results were obtained by launching, at the time of supernova 
explosions, a set of trace particles that travel passively with the 
ejecta, following them until they were in dense, potentially star- 
forming gas. 

supernova ejecta quickly become part of the star-forming 
cold medium. The lower two panels, which represent the 
same time interval as the middle panels, show exclusively 
dense (JV > 3 10 3 cm~ 3 ) and cold (typically less than 
300 K) gas, emphasizing that part of the new supernova 
yields have already been incorporated in that medium. 
This is because the kinetic energy and supernova yields 
travel together, becoming progressively admixed into the 
warm and cold phases of the ISM as the expanding su- 
pernova shells exchange momentum with the entrained 
ISM gas. Moreover, as the kinetic energy from super- 
novae is likely one of the main contributors to the su- 
personic and super- Alfvenic turbulence that causes pro- 
tostellar collapse, supernova yields are naturally part of 
the less dense medium from which pre-stellar cores form. 
Thus, the presence of freshly-synthesized 26 Al and 60 Fe 
in star-forming regions does not require external seed- 
ing such as hypothesized in models where admixing of 
short-lived radioisotopes results from the interaction of 
a supernova shock wave with an un-contaminated pre- 
stellar core (ICameron fc Truran|[l977t IBoss et al.l I2010t 
IPan et al.H2012l) . 

5. THE SOLAR SYSTEM'S INITIAL 60 FE/ 26 AL RATIO 

The level of 26 Al at the time of Solar System forma- 
tion is well-constrained from the 26 Al- 26 Mg systemat- 
ics of the oldest Solar System solids, CA Is, which define 
an in i tial 26 A1/ 27 A1 valu e of — 5 • 10" 5 (jJacobsen et all 
I2008t Larson c t al.l l20lTT) . This ratio is comparable to 
the 26 A1/ 27 A1 values of star-forming gas inferred from 
our simulations of evolving GMCs, in particular after 
10 Myr of evolution, where average 26 Al abundances are 
similar to that of the nascent Solar System (Fig. lb). 
Although there is considerable debate in the literature 
regarding the i nitial abundance of 60 Fe at the birth of the 
Solar System (iTachibana et all [200l iKrot et al.l 120091: 
iTelus et al.l I2012T ) . a recent study of pristine differenti- 
ated meteorites and chondritic components suggest that 
the initial Solar System 60 Fe/ 56 Fe may have been as low 
as -1-10" 8 (jTang fc Dauphasll2012D . At face value, this 
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Fig. 5. — Images of projected average total mass density per 
unit volume (left) and 2 ■ 10 s times the ™A1 mass density per unit 
volume (right), shortly after the explosion of a supernova. The 
box size is 40 pc. The middle two panels correspond to a time 1 
Myr later than the upper two panels, with identical normalization. 
To enhance visibility the gray scales corresponds to actual values 
raised to 0.4. The two lowermost panels show the same projections 
as the middle panels, but with all mass with number density below 
3 ■ 10 3 cm -3 removed. 

implies an initial 60 Fe/ 26 Al of ~0.002 (assuming a solar 
27 Al/ 56 Fe = 0.104. lAsplund et al.ll2009l) . which is much 
lower than the observed 7-ray flux 60 Fe/ 26 Al ratio of 
~0.2 as well as that of any star-forming gas sampled 
from our simulations (Fig. 3). We consider three possi- 
bilities for this mismatch: (i) uncertainties in the 60 Fe su- 
pernova yields, (it) contamination of the protosolar MC 
from 26 Al-rich winds of nearby massive or asymptotic 
giant branch (AGB) stars, and (Hi) the Solar System's 
initial 60 Fe/ 56 Fe value inferred from meteorites sampled 
so far is not representative of the initial Solar System's 
bulk 60 Fe abundance. 

No plausible modification of the supernovae and stel- 
lar winds yields can locally simultaneously give the ob- 
served high Solar System level of 26 A1/ 27 A1 and the 
very low level of 60 Fe / 56 Fe inferred from meteorites 
(|Tang fc D auphas 2012). This result is robust and model 
independent, in that the average ratio of the yields in the 
simulations agree with the galact ic ratio, which is well 
constrained by 7-ray observations (jWang e t al. 2007]). As 
per the lower panel in Fig. [31 local values of 60 Fe/ 26 Al be- 
low 0.10 are exceedingly rare. Thus, uncertainties in the 



supernova yields to account for the apparent low Solar 
System 60 Fe/ 26 Al ratio can be ruled out. 

Contributions to the gas subsequently forming the So- 
lar System by winds of a massive star or an AGB star 
could result in the apparent low 60 Fe/ 26 Al inferred from 
some meteoritic studies, given that 60 Fe is not efficiently 
produced in these stellar environments. However, the 
space-time window of a star-for ming region with gas of 
this low ^°Fe/ 2 6 Al ratio is small (jKastner fc Mversill994 
Williams 2010), requiring an isolated GMC pocket where 
60 Fe has decayed away, and is in the neighborhood of a 
rare, very massive star, just before the end of its brief 
lifetime, or is in the close neighborhood of a stray AGB 
star when it emits an 26 Al-rich wind. Although this can- 
not be ruled out, we estimate from our simulations the 
probability of such event to be less than 10~ 6 . 

Alternatively, the low initial Solar System abun- 
dance of 60 Fe/ 56 Fe inferred f rom differentiated mete- 
orites (jTang fc Dauphasl I2012D may not be representa- 
tive of the average initial abund ance of 60 Fe in the Solar 
System. It has been proposed (jLarsen et al.|[20Tll) that 
large-scale heterogeneity in the initial abundance of 26 Al 
may have existed throughout the inner Solar System, 
with possibly up to 80 % reduction relative to the canon- 
ical initial 26 A1/ 27 A1 ratio defined by CAIs. The 26 Al 
heterogeneity amongst inner Solar System objects has 
been asc ribed to thermal processing of molecula r cloud 
material (<Trinquier et al.ll200a iPaton et al.ll2013lh which 
resulted in preferential loss by sublimation of a more 
volatile 26 Al-rich carrier, producing residual isotopic het- 
erogeneity. Thus, CAIs represent samples of the comple- 
mentary gaseous reservoir enriched in 26 Al by thermal 
processing, which resulted in the widespread 2 Al deple- 
tions observed among the inner Solar System bodies. If 
the 60 Fe carrier phase was significantly more volatile than 
the 26 Al carrier, then thermal processing would have ef- 
fectively removed most of the 60 Fe from the inner Solar 
System solids, leading to a significant enhancement of the 

Fe concentration of the gaseous phase at that time. If 
this interpretation is correct, we predict that CAIs, pro- 
vided they are samples of a complementary gas reservoir, 
will define internal isochron relationships consistent with 
initial 60 Fe/ 56 Fe ratios of the order 10~ 6 . 

The close agreement of the Sun's abundance pat- 
tern with that of solar-twins in the M67 open cluster 
(jOnehag et al.|[20TlT ) suggest that our Solar System may 
have formed in a GMC giving rise to a similar high-mass 
cluster. The high 26 Al and 60 Fe concentrations in star- 
forming gas indicated by our simulations of a comparable 
environment suggest that the presence of short-lived ra- 
dionuclides in the early Solar System does not require 
an unusual formation scenario. Direct measurement of a 
60 Fe/ 26 Al ratio in CAIs that is comparable to that de- 
duced from our GMC simulations would give much cre- 
dence to this proposal. 
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NASA High-End Computing Capability at Ames. 
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